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Обсуждаются вопросы получения оценок координат центра, длин полуосей эллипса и их ориентации относительно системы ко- 
ординат на основании результатов измерения координат точек, расположенных на его поверхности. Показано, что процедура 
получения оценок может быть сведена к последовательности задач линейного оценивания и определения собственных векто- 
ров и собственных значений симметричной вещественной матрицы. 


В различных приложениях (см. например [1,2]) 
возникает необходимость вычисления параметров 
эллипса (длин полуосей и их ориентации относи- 
тельно выбранной системы координат) по результа- 
там определения координат точек, находящихся на 
поверхности эллипса. Эту задачу можно также рас- 
сматривать как процедуру приближения множества 
наблюденных точек эллипсом. В настоящей работе 
предложен метод решения указанной задачи. 

Каноническое уравнение эллипса размерности л, 
центр которого совпадает с точкой х 0 , имеет вид [3]: 

(х-х 0 ) г 1К(х-х 0 ) =1. (1) 

Здесь х - вектор координат точки эллипса, а 
XV матрица вида: 

\Ѵ = Ѵ Т АѴ, (2) 

где V - унитарная матрица поворота (матрица на- 
правляющих косинусов), а А — диагональная ма- 
трица: 


где л,- - длина соответствующей полуоси эллипса, 
где /=1,2,...,я. Таким образом, выражение (2) соот- 
ветствует разложению матрицы XV по собственным 
числам и собственным векторам. 

Разлагая ( 1 ) в сумму квадратичных форм, полу- 
чим: 

х г \Ѵх — 2х \Ѵх 0 + х' 0 ІѴх Г) = 1, (3) 

или 

X X %х,Х; ~ ?Х + С = 1, (4) 

/=1 у=і і=1 

т. е. уравнение, линейное относительно 
к=п(п+\)/2+п+\ неизвестных (в силу симметрич- 
ности Щ, из которых первые л(я+1)/2 суть элемен- 
ты щ <і матрицы XV, следующие л неизвестных 

Ъ і=Х м ’г х о/’ а С = х ’}Ѵх 0 . 

7=1 


Приведем выражение (4) к виду: 

X X ~ 2 Х Х А = 1 - с - (5) 

/=1 у — 1 і =1 

Если имеются координаты т ( т>к ) точек, при- 
надлежащих поверхности эллипса, то, составив для 
каждой из них уравнение вида (5), получим систе- 
му уравнений. Решая эту систему с единичным 
вектором в правой части, что равносильно умноже- 
нию обеих частей системы на 1/(1 -С), получим 
оценки коэффициентов ѵй у =ѵѵ, 7 / (1 — С7) и Ь^Ь/ІІ—С), 
из которых сформируем матрицу ІѴ=ІѴ/( I— С) и 
Ь=Ь/(\—С) вектор соответственно. Из определения 
Ъ і следует, что Ь = ІѴх 0 , поэтому, оценка х 0 может 
быть получена в виде: 

х 0 =\Ѵ 1 Ь. 

При разложении матрицы ]Ѵ на собственные 
векторы и собственные числа получим матрицу по- 
ворота Ѵи диагональную матрицу А=А/( 1-0. Вы- 
числим сумму квадратичных форм (3), подставив 
вместо матрицы В 7 полученную при решении си- 
стемы (5) матрицу ]Ѵ. 

х Т ІѴх - 2х Т \Ѵх 0 + х г 0 ІѴх 0 = 

= (х т ХѴх-2х т ХѴх 0 + х 0 (Ѵх 0 )/(\- С), 

а, учитывая (3), получим: 

х т Ш - 2х т Ш 0 + х т 0 \Ѵх 0 = 1 /(1 - С). (6) 

Вышеприведенные рассуждения обосновывают 
следующий алгоритм определения параметров эл- 
липса по координатам набора точек, принадлежа- 
щих его поверхности: 

• измерить координаты т точек, расположенных 
на поверхности эллипса; 

• составить из них систему уравнений вида (5); 

• решить полученную систему линейных уравне- 
ний с единичным вектором в качестве правой 
части; 

• составить из первых я(я+1)/2 элементов векто- 
ра решения системы матрицу Щ а из следующих 
л элементов - вектор й; 

• получить оценку вектора координат центра эл- 
липса в виде х 0 = XX ‘й; 

• вычислить матрицу собственных векторов и 
вектор собственных значений матрицы IV; 
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• получить оценку матрицы поворота эллипса от- 
носительно выбранной системы координат как 
транспонированную матрицу собственных век- 
торов матрицы IV; 

• вычислить сумму квадратичных форм 
х т Й / х-2х т Й / х 0 +х 0 г Й / х 0 , и, деля эту сумму на эл- 
ементы вектора собственных значений матри- 
цы \Ѵ определить длины полуосей эллипса, как 
корни квадратные из полученных величин. 

Для демонстрации работы алгоритма задали 

трехмерный эллипс с полуосями, равными 235, 181 
и 27. Были вычислены координаты 20 точек, слу- 
чайным образом расположенных на поверхности 
эллипса. К этим координатам был добавлен вектор 
координат центра эллипса дг 0 =(Ю, -5, 7). Затем эти 
координаты были преобразованы с использовани- 
ем матрицы поворота вида: 



0,766 

0,643 

0 

т = 

-0,583 

0,694 

0,423 


-0,272 

0,324 

0,906 


что соответствует сумме поворота системы коорди- 
нат на 25° в плоскости У2 и на 40° в плоскости ХУ. 
Исходные и преобразованные координаты точек 
представлены в табл. 1. 

Из преобразованных координат была составле- 
на система линейных уравнений вида (5), матрица 
которой приведена в табл. 2. 

Полученная в результате решения системы ме- 
тодом наименьших квадратов матрица Нравна: 


= 10 


1,235 

-1,254 

-3,337 


-1,254 

1,677 

3,977 


-3,337 

3,977 

1,144 


Матрица собственных векторов и вектор соб- 
ственных значений матрицы IV равны соответ- 
ственно: 


Ѵ = 


0,766 

0,583 


-0,272 


0,643 

-0,694 

0,324 


-3,8 -10 15 
0,423 
0,906 


и А=(1,82005-10 5 , 3,084733- ІО' 5 , 1 ,38626799- 1 0 3 ). 

Обратные значения корней квадратных из соб- 
ственных чисел равны 233,766; 180,049 и 26,858. 


Таблица 1. Координаты точек эллипса 


і Исходные 

Преобразованные | 

х, 

х 2 

х 3 

Хі 

х 2 

х 3 

-165,151 

128,758 

0,215 

-191,581 

- 21,694 

- 47,221 

- 22,412 

60,696 

25,306 

- 49,402 

30,926 

4,284 

49,491 

151,055 

- 13,745 

-36,353 

127,236 

- 69,296 

-32,692 

103,658 

21,813 

-81,356 

53,015 

- 17,039 

- 79,074 

51,832 

- 24,221 

- 74,190 

- 27,684 

-36,857 

-105,965 

15,388 

23,990 

- 86,655 

- 54,663 

22,239 

1,162 

- 44,458 

- 26,173 

43,900 

- 43,592 

2,068 

50,405 

46,042 

25,462 

14,873 

67,609 

10,618 

196,270 

-15,454 

14,669 

165,369 

115,180 

26,826 

-119,263 

58,391 

21,572 

- 121,237 

-34,137 

1,874 

164,842 

- 128,339 

-1,948 

211,571 

11,226 

59,473 

- 154,522 

-58,596 

18,369 

- 79,225 

-139,060 

48,411 

12,992 

180,438 

1,513 

- 85,576 

129,114 

- 67,885 

- 227,324 

28,733 

- 5,336 

- 179,429 

- 132,900 

- 9,979 

98,946 

123,644 

-16,111 

18,144 

139,228 

-59,856 

-100,598 

47,553 

23,347 

-101,108 

-29,090 

8,063 

-21,164 

171,643 

8,217 

- 108,437 

103,223 

-58,093 

-219,889 

- 6,678 

9,473 

- 157,128 

- 147,912 

18,408 

12,746 

-36,313 

- 26,410 

48,093 

-30,568 

-1,590 

233,026 

6,297 

-3,363 

185,754 

148,069 

1,291 


Таблица 2. Матрица системы линейных уравнений 


3.670П0 4 

470,624 

2, 230-ТО 3 

4,156-ТО 3 

9,047-ТО 3 

1,024-ТО 3 

-191,581 

-21,694 

-47,221 

2,441-ТО 3 

956,393 

18,353 

-1, 528-ТО 3 

-211,638 

132,485 

-49,402 

30,926 

4,284 

1,322-Ю 3 

1,619-10* 

4, 802-ТО 3 

-4, 625-ТО 3 

2, 519-ТО 3 

-8, 817-ТО 3 

-36,353 

127,236 

-69,296 

6.619ТО 3 

2,811 -ІО 3 

290,319 

-4, 313-10 3 

1, 386-ТО 3 

-903,301 

-81,356 

53,015 

-17,039 

5,504-ТО 3 

766,419 

1, 358-ТО 3 

2, 054-ТО 3 

2,734-ТО 3 

1, 020-ТО 3 

-74,190 

-27,684 

-36,857 

7,509-ТО 3 

2, 988-ТО 3 

494,575 

4,737-ТО 3 

-1, 927-ТО 3 

-1,216-ТО 3 

-86,655 

-54,663 

22,239 

1,927-Ю 3 

1,900-ТО 3 

4,278 

—1,914-ТО 3 

90,801 

-90,165 

43,900 

-43,592 

2,068 

221,212 

4, 571-ТО 3 

112,734 

1,006-ТО 3 

157,918 

717,845 

14,873 

67,609 

10,618 

2,735 10* 

1,327-10* 

719,615 

1,905-10* 

4,436-ТО 3 

3, 090-ТО 3 

165,369 

115,180 

26,826 

1,470-10* 

1,165-ТО 3 

3,512 

4,139-ТО 3 

-227,212 

-63,977 

-121,237 

-34,137 

1,874 

4,476-10* 

126,020 

3,537-ТО 3 

2, 375-ТО 3 

1,258-10* 

667,635 

211,571 

11,226 

59,473 

6, 277-ТО 3 

1,934-10* 

2, 344-ТО 3 

1,102-10* 

-3, 835-ТО 3 

— 6,732-10 3 

-79,225 

-139,060 

48,411 

7, 323-ТО 3 

1,667-10* 

4,608-Ю 3 

-1,105-10* 

5,809-Ю 3 

-8,765-ТО 3 

-85,576 

129,114 

-67,885 

3,219-10* 

1,766-10* 

99,586 

2,385-10* 

1,791-ТО 3 

1, 326-ТО 3 

-179,429 

-132,900 

-9,979 

329,190 

1,938-10* 

3, 583-ТО 3 

2,526-Ю 3 

-1, 086-ТО 3 

~8, 334-ТО 3 

18,144 

139,228 

-59,856 

1,022Т0 4 

846,220 

65,012 

2, 941-ТО 3 

-815,228 

-234,551 

-101,108 

-29,090 

8,063 

1,176-10* 

1,066-10* 

3,375 -ІО 3 

-1,119-10* 

6, 299-ТО 3 

— 5,997-Ю 3 

-108,437 

103,223 

-58,093 

2,469-10* 

2,188-10* 

338,854 

2,324-10* 

-2,892-ТО 3 

~2, 723-ТО 3 

-157,128 

-147,912 

18,408 

2, 313-ТО 3 

934,392 

2,527 

— 1,470-10 3 

-76,453 

48,593 

48,093 

-30,568 

-1,590 

3,450-10* 

2,192-10* 

1,666 

2,750-10* 

239,774 

191,131 

185,754 

148,069 

1,291 


Столбцы матрицы в таблице соответствуют следующим регрессорам: х, 2 , х 2 2 , х 3 2 , х, *х 2 , х, *х 3 , х 2 *х 3 , х,, х 2 , х 3 


14 






Естественные науки 


Таблица 3 . Рассчитанные координаты точек эллипса при пе- 
реносе начала координат в его центр 


X, 

х 2 

Хз 

- 191,581 

- 21,694 

- 47,221 

- 49,402 

30,926 

4,284 

- 36,353 

127,236 

- 69,296 

- 81,356 

53,015 

- 17,039 

- 74,190 

- 27,684 

- 36,857 

- 86,655 

- 54,663 

22,239 

43,900 

- 43,592 

2,068 

14,873 

67,609 

10,618 

165,369 

115,180 

26,826 

- 121,237 

- 34,137 

1,874 

211,571 

11,226 

59,473 

- 79,225 

- 139,060 

48,411 

- 85,576 

129,114 

- 67,885 

- 179,429 

- 132,900 

- 9,979 

18,144 

139,228 

- 59,856 

- 101,108 

- 29,090 

8,063 

- 108,437 

103,223 

- 58,093 

- 157,128 

- 147,912 

18,408 

48,093 

- 30,568 

- 1,590 

185,754 

148,069 

1,291 


Вычисленное по формуле (6) значение множи- 
теля 1/(1— С) равно 1 ,010702696321, а полученные с 
использованием этого множителя значения полуо- 
сей эллипса равны соответственно 235,013; 
181,010; 27,002, т. е. практически совпадают с за- 
данными. 

Таким образом, показано, что по результатам 
определения координат точек, принадлежащих его 
поверхности, можно определить параметры эллип- 
са, находящегося в общем положении. Предложен 
алгоритм получения соответствующих оценок. Ре- 
зультаты решения численного примера подтвержда- 
ют работоспособность предложенного алгоритма. 
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